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ABSTRACT 

We study spectral variability of 1 1 ultraluminous X-ray sources (ULX) using archived XMM- 
Newton and Chandra observations. We use three models to describe the observed spectra: a 
power-law, a multi-colour disc (MCD) and a combination of these two models. We find that 
7 ULXs show a correlation between the luminosity Lx and the photon index T. Furthermore, 
4 out of these 7 ULXs also show spectral pivoting in the observed energy band. We also find 
that two ULXs show an Lx — T anti-correlation. The spectra of 4 ULXs in the sample can be 
adequately fitted with a MCD model. We compare these sources to known black hole bina- 
ries (BHB) and find that they follow similar paths in their luminosity-temperature diagrams. 
Finally we show that the 'soft excess' reported for many of these ULXs at ~ 0.2 keV seems 
to roughly follow a trend L so ft oc T~ 3 - 5 when modelled with a power-law plus a 'cool' MCD 
model. This is contrary to the L oc T 4 relation that is expected from theory and what is seen 
for many accreting BHBs. The observed trend could instead arise from disc emission beamed 
by an outflowing wind around a ~ 10 solar mass black hole. 

Key words: accretion, accretion discs - black hole physics - galaxies: individual: NGC 253, 
NGC 1313, PGC 13826, PGC 23324, PGC 28757, NGC 5204, NGC 5408, NGC 6946 - 
X-rays: binaries - X-rays: galaxies. 



1 INTRODUCTION 

Ultraluminous X-ray sources (ULX) are point like, non-nuclear 
(not an active galactic nucleus, AGN) objects with X-ray lumi- 
nosities exceeding the Eddington limit for stellar mass black holes 
(StMBH), L E dd ~ 1.3 x 10 39 (M/10M Q )ergs- 1 . Many ULXs 
have been identified as accreti ng black holes because of their 
short- and long-term variability l lMifler & Colbertll2004l and ref- 
erences therein). More tha n 300 ULXs have been discovered in 
nearby galaxies (|Liull2008l) and all types of galaxies have ULXs 
jMushotzkvll200ar Although they were d iscovered with the Ein- 
stein mission almost thirty years ago (see lFabbianoi ri989, for re- 
view), the nature of these sources is not yet well understood. 

Assuming that ULXs are isotropic emitters, the most lumi- 
nous source M82 X -l has a luminosity of Lx ~ 10 41 ergs _1 
jPatruno et alj|200(3) . Therefore, for such a high luminosity ULX 
the accreting black hole should be an intermediate-mass black 
hole (IMBH) with M > 1000M Q not to exceed the Edding- 
ton limit. One argument used to support this idea is the presence 
of a ~ 0.2 keV 'soft excess' in the spectra of several luminous 
ULXs dKaaret et alj|2003t iMiller et al.ll2003h . which could be in- 
terpreted as a signature of a 'cool' accretion disc around an IMBH. 
It is also possible, however, that this excess aris es from the emis- 
sion at a large distan ce from the central source ( iKing et alj|200ll : 
IPoutanen et alj [2007?) or, maybe, is an artifact of a complicated 
absorber or emission from some thermal plasma surrounding the 
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source dGoncalves & Sorial2006l) . Another argument used in favour 
of the IMBH interpretation is the detection of low-frequency quasi- 
periodic oscillations (QP O) from ULXs in M82, Holmberg IX 
and NGC 5408 dStrohmaver & Mushotzky||2003l ; lDewangan et al.1 
l2006r . lStrohmaver et al.ll2007h . However, similar QP O frequencie s 
have been also detected from Cyg X-l dVikhlinin et al.1 [T994) , 
which is certainly not an IMBH. 

There are also several arguments against the hypothesis that 
ULXs are powered by accretion onto IMBHs. The main issue is 
the formation of such massive objects in 'normal' stellar evolution. 
Black holes of masses greater than 80M© are not expected to form 
even in the lowest metallicity enviro nments and at solar me tallic- 
ities the upper limit is about 15M© teelczvnski et alj[2009h . Fur- 
thermore, about 30 per cent of UL Xs are associated with star for- 
mation regions jMushotzky||2006l) and a small number of them 
have been associated with young stellar counterparts (Mushotzky 
2006; [Roberts et al] r2008). These facts can be naturally explained 
if ULXs are powered by accretion onto StMBHs. 

There are numerous models put forward on how the lu- 
minosity of a StMBH can exceed its Eddington limit. Truly 
sup er-Eddington models include the advective slim disc mod- 
els i Jaroszvnski et alj[l980l : lAbramowicz et al] 1 19881 : iKaw aguchi 
120031 7 which have been successfully used to describe some 
ULX spectra dVierdavanti et al. 2006). Accretion discs with 
strong outflows JShakura & Sunvaev ll973l: lKingetal.1 l200ll ; 
IKing & Puchnarewiczll2002r . |Poutanen et al j2007h can also explain 
the high luminosities of ULXs. Other plausible scenarios involv- 
ing StMBHs that can produce the high observed luminosities in- 
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Table 1. Observation log 



ULX name 
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Coordinate 
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-25:17:49 
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0110900101 0125960101 0125960201 
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0304851201 0304851101 
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— 25'20'51 
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0110Q00101 0195Q60101 0159090101 

\) 1 1 U7UU 1 \J 1 \J 1 jLjy\J\J 1\J 1 UlJZ,UiUlWl 

0104x50001 oio4.x5iooi 0104.R51901 


3931 


NGC nil 1IT X-1 


NOC1 31 3 XMM1 
1X0 7 


03:18:19.8 


-66:29:10 


3.7 


0106860101 0150980301 0150980401 
0150280601 0150281101 0205230201 

0905910101 09059104.01 0905910501 
0905910601 04.050Q0101 




lsmr 1111 ttt x.9 

, *l V J V 1 Jl J UL/A Z 


lsmr 1 1111 ymmi 
IXO 8 


03:18:22.0 


-66:36:04 


3 7 


0106X60101 0150980101 01509X0101 

\) 1 UUOUU 1U1 Ul *JVJjLO\J 1\J 1 Ul JV./.£.OLJ JU 1 

0150280401 0150280501 0150280601 

0150981101 0905910901 0905910101 

0205230401 0205230501 0205230601 
0301860101 0405090101 




IC 342 X-6 


IC 342 ULX1 
IXO 22 


03:45:55.7 


+68:04:55 


3.9 


0093640901 0206890101 0206890201 
0206890401 


2916 2917 


HoIIULX-1 


IXO 31 


08:19:29.0 +70:42:19 


3.39 


0112520601 0112520701 0112520901 














0200470101 




HoIX ULX-1 


NGC 3031 ULX2 
IXO 34 


09:57:54.1 


+69:03:47 


3.42 


0112521001 0112521101 0200980101 


3786 3787 3788 


NGC 5204 ULX- 1 


NGC5204 XMM 1 


13:29:38.6 +58:25:06 


4.3 


0142770101 0142770301 0150650301 


3933 




IXO 77 








0405690101 0405690201 0405690501 




NGC 5408 ULX-1 


NGC5408 XMM1 


14:03:19.6 


-41:23:00 


4.8 


0112290501 0112290601 0112290701 
0112291201 0302900101 




NGC 6946 X-6 


NGC 6946 ULX3 


20:35:00.7 


+60:11:31 


5.5 


0200670301 0401360301 
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4632 4633 



a Distance estimates are fromlLiu & Bregmanl fc0q5ri except for HoII ULX-1 and NGC 5408 ULX-1, where the distances are taken 
from lKarachentsev et alj 120021) and lLlu&MrrabeTl2005l) . respectively. 



elude the photon bubble model dAronsI 1 19921 : [Begelman l200lh 
and beaming of radi ation by a relativistic jet Wording etal.l f2002: 
iFreeland etail2 006). Strong beaming, however, can be ruled out in 
some cases by observations of photo-ionised nebulae surrounding 
the ULXs because the observed line properties re quire luminosities 
comparable to the observed X-ray luminosities fPakull & Mirioni 



comparable to the observed A-ray luminosities (Fakull & Mirioni 
20031: iKaaret et alJl2004 lAbolmasov et al.ll2007l: lKaaret & Corbel 



20091) . 



Because of the limited observing band of the current X-ray 
observatories capable of detecting ULXs (~ 0.2—10 keV), it is not 
surprising that many of the models mentioned above have been 
successfully used to model the observed spectra of some particu- 
lar ULX. This is mainly because many previous studies have been 
concentrate d only on modelling single observations of individual 
ULX s (e.g. lFeng & KaarefcOOalWinter et alj2006l : IStobbart et al.1 
2006). The lack of multi epoch observations can lead to misinter- 
pretation of the physical meaning of the derived spectral param- 
eters. By studying the spectral variability of ULXs, we hope to 
learn more about the physical nature of the specific spectral compo- 
nents by comparing their variations to expectations from theoretical 
models . Some previous works have studied spectral variability of 



ULXs dFabbiano et al.l2003l:lDewangan et al.l2004l : lFeng & Kaarel 
l2006d lb HRoberts et alfcOQfllBerghea et alj|2008h . but the samples 



of ULX wit h multiple observations h ave been small Interestingly, 
the study bv lFeng & Kaaretl J2006bh showed that the two ULXs in 
NGC 1313 galaxy have different types o f spectral variability. Their 
subsequent study dFeng & Kaare l l2007t) also showed that the soft 
excess in NGC 1313 ULX-2 did not follow the L oc T 4 relation 
that is expected from theory when modelled with a multi-colour 
disc (MCD). Therefore, in order to study the spectral variability of 
ULXs in better detail, we have selected a sample of 1 1 ULXs that 
have been observed multiple times with XMM-Newton and Chan- 
dra. In this paper, we present the results of spectral analysis of these 
data sets. 



2 OBSERVATIONS AND DATA REDUCTION 

We searched the data archives for high S/N XMM-Newton and 
Chandra observations of ULXs that have been observed multiple 
times. We also required that the observations were not effected by 
photon pile-up. There were 11 sources that met the requirements: 
NGC 253 X-2, NGC 253 X-4, NGC 253 X-9, IC 342 X-6, NGC 
1313 ULX-1, NGC 1313 ULX-2, Holmberg II ULX-1, Holmberg 
IX ULX-1, NGC 5204 ULX-1, NGC 5408 ULX-1 and NGC 6946 



1 Note a recent paper by Feng & Kaaret 1 2009), who consider a large sam- 
ple of ULXs. 
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X-6 (see Table 1). In the cases where s ources have mult i ple na mes, 
we followed the naming convention of iLiu & Bregmanl J2005h . 

XMM-Newton data reduction was done using the SAS ver- 
sion 7.1.0. In order to produce a homogeneously calibrated data 
set, we processed all the observations using the observation data 
files with the latest calibration files as of March 2008. The reduc- 
tion was done according to the User's guide to the XMM-Newton 
Science Analysis SystemQ Specifically, we used the event selec- 
tions "FLAG==0 && PATTERN< =4" and "#XMMEA_EM && 
PATTERN < =1 2" for the EPIC-pn and EPIC-mos instruments, re- 
spectively. Chandra observations were reduced using the CIAO 
4.0.1 software. We processed all the observations with calibration 
files from CALDB 3.4.3 following the Chandra science threads^ 

We used circular source extraction regions of varying radii 
(between 20" to 40" for XMM-Newton and 4" to 12" for Chan- 
dra). Close-by source free regions were used as a background^ 
The resulting background subtracted spectral data were binned to 
require at least 20 counts per bin to ensure adequate statistics for the 
XSPEC spectral fitting. To ensure the best possible spectral qual- 
ity, the time periods of high background were eliminated from the 
data for both XMM-Newton and Chandra observations. In the case 
of XMM-Newton the exclusion was done by removing the time pe- 
riods where the 10—12 keV full field count rate was above 0.4 cts 
s _1 and 0.35 cts s _1 for EPIC pn and mos, respectively. For Chan- 
dra, we used the ANALYZE_LTCRV.SL script on a source free event 
file to remove the high background periods. 



3 SPECTRAL ANALYSIS AND RESULTS 

We used three models to describe the observed spectra: a POW- 
ERLAW, a MCD (DISKBB model in XSPEC) and a combination 
of these two models. The effect of the interstellar medium was 
taken into account with the WABS model in XSPEC. The hydro- 
gen column density was let as a fr ee parameter, but wa s required 
to have at least the Galactic value jDickev & Lockmanlll990h . In 
some cases, none of these simple models gave statistically accept- 
able fits to the data. In most of these cases, the fits could have been 
improved by adding more model component s like gaussian s at < 1 
keV (NGC 5204 ULX-1, OBSID 3933, see lRoberts etail .20061) : 
NGC 5408 ULX-1, OBSID 0302900101) or by changing met al 
abundances of the absorber (HoII ULX-1, see lGoad etall ( 120061) ). 
We checked that our best fitting parameters did not change signifi- 
cantly by adding these components, so we decided to ignore these 
facts, because we wanted to keep our models as simple as possible. 

We fit the spectra in the 0.4—10 keV band (or up to energies 
which are not noise dominated). XMM-Newton provides spectral 
coverage down to ~ 0.2 keV, but we decided to restrict our fitting 
range to 0.4 keV because of the known calibration uncertainties be- 
low this value0 All the parameter errors in the following figures are 
given at the 90 per cent confidence level. The errors for the fluxes 
were computed using the CFLUX convolution model in XSPEC. In 
the case of the POWERLAW model, the luminosities were calcu- 
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Figure 1. The intrinsic 0.4— 10 keV band luminosity versus the photon in- 
dex obtained from the POWERLAW model fits. Different observations of a 
given object showing spectral pivoting in the observed band are connected 
by a theoretical curve {TJ. Some sources show correlation, but no spec- 
tral pivoting and some other even show anti-correlation. NGC 1313 ULX- 
1, Holmberg II ULX-1. Holmberg IX ULX-1, NGC 5204 ULX-1 show 
Lx _ r correlation with spectral pivoting in the observed band. NGC 253 
X-4, IC 342 X-6 and NGC 5408 ULX-1 show L x - T con-elation without 
spectral pivoting, and NGC 253 X-2 and NGC 1313 ULX-2 are L x - T 
anti-correlated sources. NGC 253 X-9 shows no Lx — T correlation and 
NGC 6946 X-6 is not variable enough to make firm conclusions. 



lated for the 0.4—10 keV band. The MCD component luminosities 
are bolometric. 



3.1 Power-law type ULXs and the L% — T correlation 

A simple absorbed POWERLAW model provides a good fit for most 
of the observations in our sample (see Table lAlt . We find a correla- 
tion between the luminosity and the photon index (hereafter Lx — T 
correlation) for most of the ULXs, see Fig.Q] Si milar correlations 
were also observed f or NGC 1313 ULX-1 I Feng & Kaaretll2006bl) 
and Antennae X-ll .Feng & Kaaretll2006al) . The ULXs with the 
Lx -T correlation are NGC 253 X-4, IC 342 X-6, NGC 1313 ULX- 
1, Holmberg II ULX-1, Holmberg IX ULX-1, NGC 5204 ULX-1 
and NGC 5408 ULX-1. In contrast, two ULXs of the sample, NGC 
253 X-2 and NGC 1313 ULX-2, show an anti-correlation between 



the luminosity and T (see Fig. Q] and iFeng & Kaareq|2006al lbi for 
comparison). We also find that the luminosity of NGC 253 X-9 is 
clearly not correlated with T and that NGC 6946 X-6 is not signif- 
icantly variable to make firm conclusions. 

For NGC 1313 ULX-1, Holmberg LX ULX-1 and NGC 5204 
ULX- 1 , the modelled absorption column A ^h is also correlated with 
T (see also fig. 1 in lFeng & Kaaretll2009l . for NGC 5204 ULX-1). 
If this correlation has no physical origin, then the luminosities of 
the softest spectra would be 'artificially boosted'. To evaluate this 

Guide is available at effect on the derived intrinsic luminosities, we fixed the absorption 

http://xmm.esac.esa.int/external/xmm_user_support/documentation/sas_usg/USG.pdf . , . . , „ . , . 

' u 1,, . — — ttt it-t — ti i . *. , . — n ■ ,. . — ! — n 2 -columns to their mean values, and fitted the data again. This lead to 

" Threads are available at http://cxc.harvard.edu/ciao/threads/index.html ° 

4 In the case of XMM-Newton, the background region was chosen from the somewhat worse fits for the cases where the best fitting absorption 
same EPIC-pn CCD chip as the source columns differ most from the mean, but the Lx — T correlations 

5 See EPIC status of calibration and data analysis XMM-SOC-CAL-TN- remained. Also, the observed (i.e. absorbed) luminosities L abs are 
0018. correlated with T, but the A^h — T dependence leads to a 'shallower' 
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Figure 2. All observed EF E spectra of NGC 5204 ULX-1 fitted with an ab- 
sorbed POWERLAW model. The curves show the deconvolved model spec- 
tra. Spectral pivoting occurs at ~ 5.5 keV. 



correlation between L a b s and T. Therefore, we conclude that the 
Lx — r correlations are real. 

If the Nu — r correlation is indeed artificial, then there must 
be a low-energy break in the power-law component below ~ 1 
keV in these spectra. This could be then be interpreted as a sign 
that this component arises from Comptonization of the disc photons 
in a 'corona' surrounding the source. It is also possible, however, 
that the correlation between Nu and T has a physical origin, for 
example a growing luminosity might lead to a stronger outflow and 
therefore to a denser environment. 

In general, all the sources that have luminosities permanently 
above ~ 3 x 10 39 ergs -1 show a power-law type spectra. In 
contrast, the three sources with luminosities permanently below 
~ 3 x 10 39 ergs -1 have non-power-law (or thermal) type spec- 
tra (see next subsection). An exception to this 'rule' is NGC 1313 
ULX-2 which has non-power-law type spectra at high luminosi- 
ties and power-law type spectra at lower lumino sities very much 
like th e very high state spectra of XTE J1550-564 dKubota & Pond 
|2004) . In some sources, both POWERLAW and DISKBB models give 
statistically acceptable (or equally good) fits. We assign them to 
a non-power-law type because their variability is consistent with 
Galactic black hole binaries (BHB) when modelled using DISKBB. 
Also, a simple absorbed POWERLAW model does not provide statis- 
tically acceptable fits for NGC 6946 X-6, but the DISKBB + POW- 
ERLAW model gives good fits for most observations (see Sect. 13.31 
and Table lAlt . We therefore assign this source to a power-law type. 

Four out of the 7 ULXs with the L% - T correlation (NGC 
1313 ULX-1, Holmberg II ULX-1, Holmberg IX ULX-1 and NGC 
5204 ULX-1) show also signs of spectr al pivoting. Similar pivot- 
ing has been also obser ved in Cyg X-l dZdziarski et ai]|2002l) and 



several AGNs (see e.g. Izdziarski et ail 12003 , table 1, and refer 



ences therein). The pivoting is clearly seen for NGC 5204 ULX-1 
at ~ 5.5 keV (see Fig.[2jl and the three high luminosity observations 
of Holmberg II ULX-1 seem to pivot at ~ 3.5 keV. The pivoting in 
NGC 1313 ULX-1 and Holmberg IX ULX-1 occur at the boundary 
of the XMM-Newton energy range at ~ 9 and 10 keV, respectively. 
To evaluate whether these pivots are real, we also show in Fig. Q] 
the dependences between the luminosity in a given energy band 
(.Ei, E2) and the pho ton index expected for a pivoting power-law 
( IZdziarski et alj|2003l) : 



-E 2 



AtyD z F Ei -e 2 =CE L p 



P 2-r p 2-r 

t £j 2 - Mj x 



(1) 



where D is the distance, Fe 1 -e 2 is the flux in the selected energy 
band and C is a constant. Equation l[T} gives a good representa- 
tion of the observed Lx — T relation for all these pivoting ULXs. 
We therefore conclude, that the pivoting is most likely real even in 
the cases of NGC 1313 ULX-1 and Holmberg IX ULX-1. How- 
ever, when the spectrum gets very soft at T > 3, NGC 1313 ULX- 
1, Holmberg II ULX-1 and NGC 5204 ULX-1 also show devia- 
tions from this depende nce. These deviations to the low/soft states 
jPewangan et al.l2 004). are similar to those seen in Galactic BHBs 
(see Zdziarski et al. I2002L fig. 16a). 

The process producing the observed power-law type spectrum 
is most likely Compton up-scattering of soft disc photons in a hot 
plasma surrounding the inner part of the disc. The observed Lx — T 
correlation and the pivoting in the sources can be explained by vari- 
ability of the luminosity of the soft disc p hotons. A larger disc l urni- 
nosity leads to a softer X-ray spectrum (Zdziarski et al. 2003|). An 
example of such variability for a cons t ant ho t-plasma luminosity is 
given in fig. 3 of IZdziarski & Grandil d200ll) for 3C 120. The sim- 
ilarity of this source to NGC 5204 ULX-1 (and the other 3 ULXs 
with pivoting) is remarkable (see Fig. [2j . Transition to the low/soft 
state can be associated with the decreasing luminosity of the hot 
plasma (corona) with the nearly constant disc luminosity. 

The highest quality data also seems to s how a break/cutoff in 
the power-law compo nent above ~ 3 keV dStobbartet ail [2006; 
iGladstone et~alT 2009). This coul d be then interpreted a s a sign of 
a cool and optically thick corona l IGladstone et al ]|2009h . Unfortu- 
nately the limited observing band of XMM-Newton and Chandra 
and the typical data quality in our sample do not allow us to see 
these high energy breaks/cutoffs in most of these data. A dedicated 
monitoring campaing of these power-law type ULXs would help 
us to check whether these breaks are always present, which would 
be highly useful in order to understand the nature of the emission 
mechanisms in ULXs. 



3.2 Non-power-law type ULXs 

Out of the selected 11 ULXs, four ULXs can be adequately fit- 
ted with an absorbed MCD model. These sources include the two 
ULXs that show the L x - T anti-correlation (NGC 253 X-2 and 
NGC 1313 X-2) and the two low-luminosity ULXs (NGC 253 X-4 
and NGC 253 X-9). The resulting fits are shown in a luminosity- 
temperature (hereafter L — T) diagram in Fig. [3] together with 
the ob servations on known BHBs adopted from lGierliriski & Pond 
d2004l) . These four ULXs follow the behaviour seen in BHBs. The 
properties of the sources can be summarised as follows. 

(i) NGC 253 X-2 closely follows the L cx T 4 relation. Its track 
on the diagram is a high luminosity-high temperature extension to 
XTE J1550-564 and LMC X-3. The relatively high luminosities 
and hot inner disc temperatures of this source compared to BHBs 
could be explained by assuming that this source is a maximally 
rotating Kerr black hole with M ~ 70M© and inclination of ~ 70 
degrees dHui & Krolikll2008l) . 

(ii) NGC 253 X-4 on the other hand behaves like LMC X-l which 
cl early does not fo l low th e expected L cx T 4 relation. As discussed 
in lPoutanen et al.l d2007h . this can be understood if these sources 
are StMBHs (or even neutron stars) accreting slightly above their 
Eddington accretion rates. In this case the observed temperature 
corresponds to the emission from the spherization radius or even 
the distant photosphere if observed edge-on. 

(iii) NGC 253 X-9 first seems to follow the L cx T 4 track at lower 
luminosities, but then at L > 10 39 ergs -1 deviates towards lower 
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Figure 3. The MCD model luminosity-temperature diagram for the non- 
power-law type ULXs. All the sources roughly follow the tracks of known 
BHBs shown by small symbols. 



temperatures similarly to LMC X-l and NGC 253 X-4. This type 
of behaviour is also seen in the spectra of XTE J1550-564 when 
it is in the very high state (VHS) instead of the high/soft state 
dKubota & Donell2004l) . 

(iv) NGC 1313 X-2 shows a slight deviation from the L oc T 4 
track much like XTE J1650-500. As the luminosity of this source 
greatly exceeds the Eddington limit for a 10 solar mass black hole, 
the simple DISKBB model is not anymore a good approximation of 
the spectrum. Indeed the DISKBB model failed to properly model 
the high-energy curvature of the spectrum, predicting slightly too 
steep cutoff. Because at lower luminosities the observed spectra 
of this source is very much like VHS spectra of XTE J1550-564 
dKubota & Donell2004l) . the higher luminosity non-power- law type 
spect r a could be understood as an ultraluminous state dRobertsI 
120071: ISoriall2007l) . However, a more elaborate model than the 
DISKBB should be used to investigate this state further. 



3.3 The soft excess and the nature of ULXs 

When fitting the data with the POWERLAW plus a cool MCD model, 
the fit improves for many (but not all) observations compared to the 
single component fits (see Table lAlV The L so ft — Tin relation for 
this soft excess is shown in Fig. [4] We find that the intrinsic lu- 
minosity of the MCD component roughly scales as L BO f t oc T-~ 3 ' 5 . 
The b est fitting relation derived using the fitexy routine dPress et all 
1 19921) is L soft oc T;~ 3 5±0 - 2 with xLi ~ 1-09. However, this result 
should be taken with some caution and there are several points to 
be considered. 

(i) The highly luminous low-temperature excesses (L so ft > 
10 40 ergs" 1 , Ti n ~ 0.1 keV) of NGC 253 X-2 do not seem to 
be real features, but rather artifacts of improper modelling of the 
harder emission with a power-law. This results in a rise in the mod- 
elled absorption column, which then causes the increase in the esti- 
mated intrinsic luminosity of the soft excess (see Fig. [5}- We there- 
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Figure 4. The luminosity-temperature diagram of the soft excess from the 
POWERLAW plus cool MCD model fits (only one set of errors per source 
is shown for clarity). All the sources with this model follow roughly a 
L so ft oc Tj~ 3,5 scaling. Note however that this scaling is subject to some 
modelling uncertainties (see subsection [33) and that the coolest and high- 
est luminosity soft excesses of NGC 253 X-2 are artifacts of an improper 
modelling of the hard emission (see Fig.|5]and Table[AT). 
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Figure 5. EF E spectra of NGC 253 X-2 fitted with an absorbed POWER- 
LAW plus cool MCD model. EPIC pn, mosl and mos2 data are shown in 
blue, red and green, respectively. The MCD component is shown by dot- 
dashed line. The dashed line shows this with the effect of absorption re- 
moved. A simple absorbed MCD fit with T ln ~ 1.5 keV improves the fit 
from Xr e d ~ 1.17 to X^ed ~ 1-08- Improper modelling of the continuum 
by a power-law leads to almost an order of magnitude increase in the ab- 
sorption column density, which produces the highly luminous cool, ~ 0.1 
keV, MCD component. 



fore did not include the soft excesses of NGC 253 X-2 in deriving 
the best fitting L so f t — Ti n relation. 

(ii) As in the POWERLAW fits, 7V H and V are still correlated for 
NGC 1313 ULX-1 and NGC 5204 ULX-1, but no longer for Holm- 
berg IX ULX-1. As discussed earlier (see subsection I3.lt . there 
is no clear evidence whether this correlation has a physical ori- 
gin, or the power-law component has a low-energy break, what one 
would expect at about Tin, if it arises from Comptonization of the 
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disc photons in a corona. If such break would exist in the power- 
law component (note that Y > 2 for NGC 1313 ULX-1 and NGC 
5204 ULX-1), then it would necessarily lead to both lower absorp- 
tion columns and higher absorbed MCD component luminosities to 
compensate the 'loss' of flux at low energies. But on the other hand, 
the decrease in TVh could then ultimately lead to lower deabsorbed 
luminosities L so ft • 

(iii) If one interprets the spectra in terms of an outflow model 
i Poutanen et al .1120071) . one would then expect an increase of TVh 
with the luminosity because the outflow would increase the column 
depth. Also, Nn could be then inversely proportional to Xl n , as it 
should decrease with increasing luminosity. 

(iv) The host galaxy distances have some uncertainties, which can 
lead to systematic errors in the derived luminosities. 

(v) The MCD component appears at low energies, where the cali- 
bration of XMM-Newton and Chandra are the least accurate. Also, 
at these energies the correct estimation of the absorption column 
makes the most difference in the derived parameter values of the 
MCD component. 

(vi) In addition to these complications, the selection of the lower 
energy bound in the spectral fitting also has some effect. For exam- 
ple, there are some dis crepancies be t ween our results of NGC 5204 
ULX-1 and those of lFeng & Kaaretl d2009l) . who used 0.2 keV as a 
lower energy bound in their fitting. 

All these factors cast a doubt on the derived best fitting relation 
L so ft tx y.- 3 ' 5 ± - 2 > 5 u t the overall trend of cooler discs with higher 
luminosities at L so ft < 10 40 erg s _1 and T ln > 0.1 keV seems to 
be real. However, the observed trend (whether or not it has a power- 
law index —3.5 or something close) is not consistent with the in- 
terpretation of an accretion disc around an IMBH. If the observed 
soft excess does not originate from a disc around an IMBH, what is 
its origin? One possible explanation could be that the thermal emis- 
sion arises from t he spherization radius of a super-Eddington disc 
with an outflow dPoutanen et alj|2007h . As observed, this model 
predicts that the thermal emission becomes cooler with increasing 
luminosity. Although the o bserved trend is too s teep for the pre- 
dicted relation of this model ( Poutanen et al .120071 . eqs. 19, 36-38), 
the apparent luminosity of this component might be overestimated 
because the out flowing wind causes some of the emission to be 
beamed. Indeed, King & Puchnar ewicj d2002l) showed that at low 
inclination angles an outflow causes the observed luminosity to be 

^-2.3xl0-(^)" 4 ^ergs-, (2) 

where the luminosity is scaled to the Eddington value / = L/L-^dd, 
p is a factor allowing for deviation from the spherical symmetry, 
r is the radius scaled to the Schwarzschild radius r — R/R s 
(Rs = 2GM/ c 2 ) and b is the beaming factor. Assuming that 
Lsoft — Lobs, King (2009) argued that the beaming factor has 
to scale as b tx m~ 2 , where the accretion rate is scaled to its 
Eddington value m = M/A/ecM, to satisfy t he observed relation 

-^soft 

~ 7 x 10 40 T ; ; 4 based on the res ults of iKaiava & PoutanerJ 
(2008). Using this scaling, iKingl d2009h estimated the mass of the 
accreting black hole 

M 4500 L soft 

— x (3) 

M m 2 (l + lnm) lO^ergs- 1 ' y 

where the factor x is the order of unity (see lKindl2009l. eq. 7). 
In our view, however, the estimation of the black hole mass given 
by equation l[3} must be taken with caution, because of the possi- 
ble systematic uncertainties in the derived best fitting parameters 



of the MCD component. Never the less, we conclude that our ob- 
servational data is consistent with the idea that the most luminous 
sources in our sample are ~ IOM0 black holes accreting at mildly 
super-Eddington rates m ~ 10 and viewed at low inclinations. 

The ULXs that have this soft excess in our sample have a 
power-law type spectra. This spectrum most likely arises from 
Comptonization of the disc photons in a corona surrounding the 
central parts (within the spherization radius) of the acc retion flow 
dHaardt & Maraschil \l99J : IPoutanen & Svenssonlll996l) . Another 
possibility for the power-law spectra could be due to overheating o f 
the accreting gas close to accreting black hole dBeloborodovll 19981) . 
Therefore, a plausible explanation to both the observed L so ft — Tin 
relation and power-law type spectra could be that at low inclina- 
tions the central Comptonized part of the disc would be visible 
to the observer at all times and the outflowing wind would cause 
the Comptonized emission and the emission from the spheriza- 
tion radius to be geometrically beamed, thus resulting in a high 
~ 10 40 ergs -1 apparent luminosity. On the other hand, at higher 
inclinations the central part of the disc would be obscured by the 
outflow. This would result in a non-power-law type spectra that we 
see for NGC 253 X-4 and NGC 253 X-9. 



4 CONCLUSIONS 

We have studied spectral variability of 1 1 ULXs and can draw sev- 
eral conclusions. 

(i) The ULXs of this sample can be roughly separated into two 
types based on their spectra: power-law- and non-power-law type 
ULXs. All the power-law type ULXs are constantly more lumi- 
nous than ~ 3 x 10 39 ergs -1 , i.e. they show permanently super- 
Eddington luminosities for a ~ 10M© black hole. 

(ii) Six out of the seven power-law type ULXs show a Lx — T 
correlation, i.e. the observed spectrum softens as the luminosity in- 
creases. Furthermore, four ULXs show also spectral pivoting with 
pivot energies E p < 10 keV. Three ULXs (NGC 1313 ULX-1, 
Holmberg II ULX-1 and NGC 5204 ULX-1) also show state tran- 
sitions to a low/soft state. 

(iii) The non-power-law type ULXs can be adequately modelled 
with a MCD. These sources also follow the behaviour of known 
BHBs in the L — T diagrams. 

(iv) NGC 1313 ULX-2 is a very interesting source showing a VHS 
type spectra (similar to XTE J 1550-564) at lower luminosities and 
a non-power-law (thermal) like spectra at high luminosities. It is the 
only source in this sample showing transitions between a spectral 
state that is seen in known BHBs and an ultraluminous spectral 
state. NGC 1313 ULX-2 is therefore most likely a super-Eddington 
source. 

(v) The soft excess seems to roughly follow a trend L so ft oc 
T ; ~ 3,5 , which does not support the idea that the observed soft ex- 
cess arises from the disc emission around IMBHs. Instead, the soft 
excess in consistent with beamed disc emission from the spheriza- 
tion radius. 

We therefore conclude that the ULXs are most likely pow- 
ered by super-Eddington accretion onto ~ IOMq black holes. The 
two different types of ULXs result from their different inclination 
angles. The permanently super-Eddington power-law type ULXs 
might be viewed at low inclination angles. The low-luminosity 
ULXs with non-power-law spectra are consistent of being super- 
Eddington sources viewed at high inclinations (NGC 253 X-4, 
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NGC 253 X-9) or they are sub-Eddington sources (NGC 253 X- 
2). 
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APPENDIX A: BEST-FITTING MODEL PARAMETERS 



Table Al: Best-fitting parameters and absorption corrected luminosity estimates. 
The results are given only if x? e d < 2. In the case of the POWERLAW + DISKBB 
model, the results are given if the soft excess DISKBB component is constrained. 
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a Hydrogen column density in units 10 21 cm 2 . 
b Photon index. 

c Inner disk temperature of the DISKBB component in units of keV. 

d Power-law model luminosity in 0.4-10 keV band in units of 10 39 ergs -1 . 

e DISKBB model bolometric luminosity (computed in 0.01-50 keV band) in units of 10 39 ergs -1 . 

f Reduced \ 2 and the degrees of freedom. 

9 Rejection probability. 



